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ABSTRACT 

There  are  some  indications  that  in  large  measure  the  density  anomaly  field  of  the 
California  Current  System  (CCS)  follows  temperature,  salinity  playing  a  significantly 
lesser  role.  These  indications  have  been  used  as  justification  for  ignoring  salinity  vari- 
ations in  dynamic  models.  An  extensive  data  base  of  simultaneous  temperature-salinity 
observations  taken  in  the  CCS  is  used  to  calculate  and  assess  the  saline  contribution  to 
specific  volume  anomaly,  dynamic  height,  and  density  variance.  The  results  show  that 
the  distribution  of  salinity  can  be  important  in  defining  the  large  scale  circulation  of  the 
CCS,  and  that  the  local  variability  can  be  quite  significant. 
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I.     INTRODUCTION:   IMPORTANCE  OF  SALINITY 

At  a  given  pressure,  the  density  of  seawater  depends  on  temperature  and  salinity. 
Except  for  high  latitude  regions  and  near  sources  and  sinks  of  salt,  temperature  is  usually 
considered  to  van.-  more  than  salinity  in  the  upper  ocean,  and  is  generally  considered  to 
be  the  major  source  of  density  variations  (Pickard  and  Emery,  1982).  Recently  Cooper 
(19S8)  has  demonstrated  the  importance  of  salinity  as  a  consideration  in  modeling 
motions  in  the  tropics.  This  result  is  surprising  because  in  the  tropics  the  thermal  ex- 
pansion coefficient,  which  increases  with  temperature  and  salinity,  is  large,  and  the 
salinity  contraction  coefficient,  which  decreases  as  temperature  and  salinity  increase,  is 
relatively  small,  and  temperature  variation  might  reasonably  be  expected  to  dominate 
forcing  of  density  differentials.  Is  it  possible  that  salinity  variability  in  the  California 
Current  System  (CCS),  which  is  geographically  between  the  polar  regions  and  the  trop- 
ics, is  also  an  important  consideration  for  modeling  dynamics?  A  close  examination  of 
the  southern  CCS  area  in  Figure  1  suggests  that  salinity  variation  could  have  as  im- 
portant an  effect  in  the  CCS  as  in  the  tropics. 

The  CCS  is  in  the  eastern  quadrant  of  the  North  Pacific  anticyclonic  gyre,  and,  ac- 
cording to  Hickey  (1979),  is  composed  of  four  currents.  The  California  Current  (CC)  is 
on  the  surface  and  extends  as  far  as  1000  km  from  the  coast.  A  surface  poleward  flow 
exists  during  the  fall  and  winter  north  of  Point  Conception  called  the  Davidson  Current. 
South  of  Point  Conception,  and  inshore  of  the  Channel  Islands  is  another  surface  How 
called  the  Southern  California  Countercurrent.  Finally,  the  California  Undercurrent 
(CU)  is  a  subsurface  poleward  flow  that  occurs  over  the  continental  slope.  Typical  mean 
velocities  associated  with  the  CC  are  on  the  order  of  10  cm  srx,  but  superimposed  on  the 
mean  flow  are  small  scale  eddies  and  jets  (Bernstein  et  al,  1977;  Mooers  and  Robinson, 
1984;  Rienecker  et  al,  1985,  1988).  The  meanders  can  have  wavelengths  of  several 
hundred  kilometers,  can  intensify  and  separate  from  the  mean  flow  and  can  exist  as 
isolated  cyclonic  or  anticyclonic  eddies  (Bernstein  et  al,  1977).  Baroclinic  jets,  typically 
~70  km  wide,  with  peak  velocities  near  80  cm  s~l  are  embedded  in  the  eddy  field  ( Mooers 
and  Robinson,  1984;  Rienecker  et  al.,  1985;  Flament  et  al,  1985;  Kosro  and  Huyer, 
1986).  Mean  surface  temperatures  in  the  CCS  range  roughly  from  10  to  20°  C  while 
surface  salinities  increase  from  32  psu  in  the  north  to  34  psu  in  the  south. 
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Figure  1.  Dynamic  Height  of  (he  Tropical  Ocean  0/400  in:  Calculated  from 
Lcvitus  (1982)  data,  (a)  Salinity  set  to  constant  }S,  (b)  salinity  varying 
as  in  Lcvitus  data,  and  (c)  difference  between  (a)  and  (b).  Contour  in- 
tervals arc  10  dyncm  in  (a)  and  (b),  and  5  dyncm  in  (c)  (From  Cooper. 
1988). 


A  typical  section  across  the  CCS.  from  San  Francisco  Bay  southwest  to  33°  X, 
134°  W  (Figure  2).  shows  that  isopycnals  very  closely  follow  isotherms,  and  intersect 
isohalines.  The  implication  is  that  temperature  variation  is  much  more  important  than 
salinity  variation  in  driving  the  density  gradient.  On  the  other  hand,  examination  of 
Figure  3,  which  depicts  the  coastal  area  of  the  CCS  just  north  of  Point  Arena  (at 
~39°  N  and  near  the  shoreward  end  of  Figure  2)  shows  that  at  a  depth  of  100  m,  about 
a  third  of  the  cross-shore  density  gradient  over  a  degree  of  longitude  is  due  to  salinity 
change.  In  the  CCS,  unlike  all  other  eastern  boundary  currents,  salinity  increases  with 
depth  (Wooster  and  Reid,  1963).  Therefore,  if  this  cold,  salty  100  m  water  were  upwelled 
to  the  surface,  the  resulting  horizontal  density  gradient  would  be  one  third  greater  than 
if  the  salinity  profile  were  constant. 

The  approximate  equation  of  state  for  seawater  is: 

P  =  P0ll-ct(T-T0)  +  (S{S-S0)l  (1) 

which  can  be  rewritten  as: 

Ap  =  -ap0Ar+/?p0AS,  (2) 

where  a  and  /?  are  the  expansion  and  contraction  coefficients  for  temperature.  T.  and 
salinity,  S,  respectively,  p  is  density,  S0  is  the  reference  salinity  (35  psu)  and  Tn  is  the 
reference  temperature  (0°  C). 

This  equation  can  be  employed  to  get  a  rough  idea  of  the  large  scale  importance 
of  salinity  if  typical  CCS  values  for  T  and  S  are  used.  Entering  the  International 
Oceanographic  Tables  (UNESCO,  1987)  for  a  and  (3  with  15°C  and  33  psu,  and  con- 
sidering typical  surface  temperature  and  salinity  ranges  of  10°  C  and  2  psu  for  the  CCS, 
computations  of  the  contributions  of  temperature  and  salinity  to  the  density  anomaly 
are  10.5  and  7.5  kg  nr3,  respectively.  Annually,  a  mean  temperature  of  12°  C  plus  or 
minus  3°  C  variation,  associated  with  a  mean  salinity  of  33  with  .3  psu  variation,  is 
roughly  characteristic  of  the  central  region  of  the  CCS.  Use  of  the  1987  tables  here 
shows  a  1.7  and  2.7  kg  nr3  contribution  to  the  annual  density  anomaly  change  for 
salinity  and  temperature,  respectively,  over  the  course  of  a  year.  Vertically  a  mean 
temperature  at  500  m  of  about  5°  C  plus  or  minus  5°  C  surface  to  1000  m,  in  conjunction 
with  an  approximate  mean  salinity  of  34,  plus  or  minus  0.5  psu,  yields  density  anomaly 
changes  of  7.2  kg  irr3  from  temperature  and  3.8  kg  m-3  from  salinity. 


These  crude  calculations  seem  to  indicate  that  the  importance  of  salinity  variation 
should  not  be  summarily  dismissed  when  considering  CCS  density  gradients  in  any  di- 
mension. Yet  often  dynamic  models  of  the  CCS  (e.g.,  Batteen  el  al.,  1989)  do  not  allow 
salinity  contributions  to  density  gradients.   What  is  the  price  of  this  omission? 

The  objective  of  this  thesis  is  to  address  that  question  by  investigating  the  saline 
impact  on  the  CCS.  The  data  set  used  is  described  in  Section  II.  A  general  description 
of  the  characteristic  salinity  of  the  CCS  is  provided  in  Section  III.  In  Section  IV  several 
techniques  are  used  to  analyze  the  data.  Mean  values  for  each  of  the  constituents  of  the 
thermosteric  anomaly,  SrSs  and  Ssn  are  computed  at  various  depths  and  analyzed.  Mean 
dynamic  heights  are  computed  from  salinity  and  temperature  observations  and  are 
compared  to  dynamic  heights  computed  from  temperature  and  constant  salinity  values. 
The  equation  of  state  for  seawater  is  differentiated  with  respect  to  temperature  and 
salinity,  variances  and  covariances  of  temperature  and  salinity  are  computed,  and  the 
contributions  of  each  to  variance  in  the  density  field  are  determined.  A  summary  of  the 
results  is  presented  in  Section  V  and  recommendations  for  future  work  are  in  Section 
VI. 
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Figure  2.  CalCOFI  Line  60,  July:  33°  N,  134°W,  to  37"  N,  123°  W:  (a)  Tem- 
perature (°C),  (b)  Salinity  (ppt),  (c)  Density  Anomaly  (from  Lynn  el  a!., 
1982). 
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Figure  3.  Temperature,  Salinity,  and  Density  Gradients:  Near  Pt.  Arena  (39°  N), 
at  100  m  depth:  (a)  Temperature  (°  C),  (b)  Salinity  (ppt),  (c)  Density 
Anomaly  (from  Iluycr  and  Kosro,  1987). 


II.     DATA  DESCRIPTION 
A.     SOURCES  AND  TYPES  OF  DATA 

The  data  base  used  for  all  calculations  was  compiled  by  Dynalysis  of  Princeton,  and 
is  described  in  detail  in  their  final  report  to  the  U.S.  Department  of  the  Interior 
(Blumberg  el  ai,  1984).  Coverage  extends  along  the  west  coast  of  North  America  to 
133  °  W.  and  from  20°  to  50°  N.  Sources  include  the  National  Oceanographic  Data 
Center  (NODC)  historical  data  archives,  the  Fleet  Numerical  Oceanography  Center 
(FNOC)  data  files,  California  Cooperative  Oceanic  Fisheries  Investigation  (CalCOFI) 
surveys  which  have  been  made  systematically  since  1949,  and  the  Coastal  Ocean  Dy- 
namics Experiment  (CODE).  CODE  data  are  the  most  current,  dating  from  as  recently 
as  1983,  and  cover  the  coastal  area  near  Point  Arena  within  a  one  degree  square  domain. 
The  highest  density  of  NODC  data  is  the  result  of  observations  from  Oregon  State 
University  standard  surveys  along  the  "Newport"  line  (~  45°  N)  and  University  of 
Washington  surveys  of  the  Columbia  River  plume  (46-47°  N).  Original  files  overlapped. 
but  the  data  have  been  carefully  screened  by  Dynalysis  to  eliminate  duplicate  observa- 
tions. The  combined  data  set  consists  of  over  300,000  data  points  with  observations  at 
up  to  33  standard  levels  to  4000  m.  The  more  than  40,000  of  these  points  which  included 
any  simultaneous  temperature  and  salinity  observations  at  500  m  or  above,  were  con- 
sidered here,  while  those  including  only  temperature  were  not  used.  Original  observa- 
tions of  salinity  were  reported  in  units  of  parts  per  thousand  (ppt)  rather  than  practical 
salinity  units  (psu).  However,  because  the  difference  in  units  is  very  slight  (Lewis  and 
Perkin,  1981)  and  not  significant  at  the  scales  discussed  here,  the  currently  preferred  psu 
is  used  throughout  this  thesis.  Most  of  the  retained  records  are  from  CalCOFI  surveys. 
which  were  collected  at  average  intervals  of  74  km  along  lines  spaced  74  km  apart  cov- 
ering from  San  Francisco  to  Baja  California  (Reid,  1988).  As  a  result  the  CalCOFI 
collection  defines  for  the  most  part  the  maximum  resolution  possible  from  the  data  set, 
exceptions  being  along  the  Oregon  coast,  and  near  Point  Arena,  where  finer  resolution 
may  be  possible.  Generally  the  number  of  observations  increases  toward  the  coast,  with 
greatest  concentration  along  the  coasts  of  California  and  Oregon.  There  are  very  few 
data  more  than  500  km  off  the  coast  or  between  38°  and  42°  N.  Discussion  here  is 
limited  to  the  area  bounded  by  the  extent  of  this  data. 


B.     DATA  MANIPULATION 

For  this  study,  data  were  placed  in  bins  identified  by  one  degree  squares  centered 
at  half  degree  latitude-longitude  intersections,  by  month  or  season,  and  by  depth  level. 
Information  so  cataloged  is  referred  to  as  a  data  block  throughout  this  paper.  Standard 
depth  levels  0,  10,  20.  30,  50,  75,  100,  125,  150,  200.  250,  300,  400,  and  500  m  were 
considered.  Choice  of  500  m  as  the  maximum  depth  level  was  made  due  to  data  avail- 
ability and  because  the  phenomenon  of  interest,  the  CCS,  has  been  classically  considered 
shoaler  than  that  depth,  at  least  within  300  km  of  shore  (Hickey,  1979).  Choice  of  one 
degree  squares  was  made  again  due  to  data  availability;  the  CalCOFI  sampling  grid  does 
not  support  finer  resolution. 

Plots  of  the  numbers  of  observations  associated  with  each  data  block  for  the  months 
January,  April,  July  and  October  are  in  Appendix  A.  When  mean  values  and  variances 
were  calculated,  at  least  ten  observations  in  a  data  block  were  required,  except  in  the 
calculation  of  mean  dynamic  height,  where  five  observations  were  used.  Mean  values 
generated  from  these  numbers  of  observations  were  considered  statistically  significant 
and  used  unweighted  in  subsequent  analyses.  Weighting  based  on  redundancy  was  not 
applied  to  avoid  biasing  results  in  favor  of  areas  of  frequent  observation.  The  determi- 
nation often  observations  as  statistically  significant  was  done  subjectively  after  studying 
the  data.  Unfortunately  specific  years  over  which  statistics  for  a  given  data  block  were 
not  considered,  and  a  bias  toward  a  particular  year,  or  years,  may  well  exist.  For  mean 
dynamic  height  calculations,  only  casts  which  include  temperature  and  salinity  data  at 
all  standard  depths  from  the  surface  to  500  m  were  used. 


III.     CHARACTERISTIC  SALINITY  OF  THE  CALIFORNIA  CURRENT 
SYSTEM 

A.  GENERAL  FEATURES 

Along  the  west  coast  of  the  United  States,  sea  surface  salinity  generally  increases 
from  about  32  psu  near  the  Washington-Canadian  border  to  about  34  psu  near  Baja 
California.  At  a  given  latitude  salinity  tends  to  be  maximum  at  the  coast,  decrease  off- 
shore, then  increase  again  at  the  western  limits  of  the  CCS  (Figure  4).  That  is,  the 
isohalines  describe  a  trough,  the  axis  of  which  extends  from  the  coast  near  the  mouth 
of  the  Columbia  River  near  46°  15'  N  to  a  point  about  500  km  off  the  coast  at  Baja 
California.  In  the  southwest  corner  of  the  CCS,  approximately  700  km  off  the  coast,  the 
isohalines  are  oriented  northwest-southeast  with  a  relatively  tight  gradient,  of  the  order 
0. 1  psu  in  40  km. 

As  mentioned  earlier,  the  stratification  of  the  CCS.  unlike  that  of  other  eastern 
boundary  currents,  has  salinity  generally  increasing  with  depth.  Profiles  (Figure  5)  are 
marked  by  permanent  haloclines:  sharp  and  shallow,  150  m  (Huyer,  1983)  in  the  north 
and  less  distinct  and  deeper,  200  m.  to  the  south  (Emery  and  Dewar,  1982).  At  500  m 
the  mean  salinity  increases  from  34.1  psu  near  50°  N  to  34.5  psu  near  20°  N. 

B.  VARIABILITY 

The  most  dramatic  seasonal  change  to  the  surface  salinity  in  the  CCS  results  from 
river  discharge.  During  the  summer  the  only  important  sources  of  fresh  water  runoff 
along  the  West  Coast  are  the  Fraser  River,  which  discharges  through  the  Straits  of 
Georgia  and  Juan  de  Fuca  at  48°  30  '  N,  and  the  Columbia  River  at  46°  15v  N.  The 
Fraser's  fresh  water  is  mixed  deep  and  its  signature  lost,  but  the  direct  effect  of  the 
Columbia  is  visible  as  far  south  as  San  Francisco  (Huyer,  1983).  During  spring  an  area 
of  tight  isohaline  spacing,  roughly  0.1  psu  in  30  km,  exists  within  ~500  km  radius  of  the 
Columbia's  mouth.  Through  the  course  of  a  year  the  area  of  tight  gradient  (Figure  4) 
gradually  decreases  until  by  winter  the  strong  gradient  exists  only  in  the  immediate  vi- 
cinity of  the  mouth  of  the  river  and  northward  along  the  coast  of  Washington. 

Using  CalCOFI  data  collected  as  far  north  as  San  Francisco  through  1978,  Lynn 
and  Simpson  (1987)  have  examined  the  physical  characteristics  of  the  CCS  including 
salinity.  The  salinity  minimum  is  at  the  surface  off  the  coast  of  northern  California  and 
50  to  100  m  deep  in  the  Baja  California  region.     Near  southern  California  there  is 
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Figure  4.      Surface  Salinity  of  the  CCS:     Units  are  ppt.  (a)  Summer  (b)  Winter 
(from  Blumberg   et  a/.,  1984). 


Salinity  (psu) 

32 

33 

33 

34 

0   m          1 

r 

100  m 

\ 

300  in 

\ 

500  in 

46°  N  t  127°  W. 

26°  /vIf  ur  w. 

FEBRUARY 

Figure  5.  Typical  Salinity  Profiles  of  the  CCS:  These  salinity  (psu)  profiles  dis- 
play mean  values  for  February  generated  from  the  CCS  data  base  used 
for  this  study. 


generally  a  minimum  at  the  surface,  but  in  the  autumn  a  subsurface  minimum  appears. 
Lynn  and  Simpson  (1987)  have  noted  significant  correlations  between  the  location  of  the 
salinity  minimum  and  the  behavior  of  the  CC.  In  January,  a  strong  inshore  poleward 
current  exists  while  the  salinity  minimum  near  Baja  California  is  below  the  surface  and 
offshore  of  the  continental  slope.  As  this  minimum  moves  inshore  and  to  the  surface, 
the  inshore  flow  becomes  strongly  equatorward.  In  October,  the  salinity  distribution  is 
similar  to  that  of  January,  but  the  inshore  flow  is  equatorward.  Elsewhere  in  the 
CalCOFI  region,  the  area  of  minimum  salinity  is  consistently  near,  but  seaward  of.  the 
maximum  equatorward  flow  of  4  to  20  cm  s~l.  Because  of  the  high  correlation  of  salinity 
minimum  to  velocity  maximum  to  a  depth  of  100  m,  Lynn  and  Simpson  (19S7)  suggest 
that  the  low  salinity  core  can  describe  the  path  of  the  CC. 


IV.     EFFECT  OF  SALINITY 

A.  STRATEGY 

Several  techniques  were  used  to  analyze  the  data  and  are  described  in  the  following 
sections:  specific  volume  anomaly  (Section  IVB),  variance  and  correlation  of  the  specific 
volume  anomaly  constituents  (Section  IVC),  mean  differences  of  the  specific  volume 
anomaly  constituents  (Section  IVD),  dynamic  height  (Section  IVE),  and  propagation 
of  the  temperature-salinity  variance  (Section  IVF).  Each  method  was  applied  at  depths 
10,  100,  200,  and  500  m.  Except  in  Section  IVF,  twelve  mean  months  were  examined. 
In  Section  IVF  the  statistical  techniques  invoked  called  for  use  of  maximum  number  of 
observations,  so  the  data  were  partitioned  into  three-month  seasons;  winter,  for  exam- 
ple, was  defined  as  January-March.  Graphics,  except  where  otherwise  noted,  were  gen- 
erated for  all  four  depths,  each  time  period,  and  each  technique.  Practical  limits  on 
space  prevent  inclusion  of  all  these  plots  in  this  thesis,  l  Only  representative  examples, 
which  when  contoured  were  hand  contoured,  are  reproduced  in  the  body  of  the  text.  In 
particular,  the  month  of  July  is  often  chosen  as  the  illustrative  period  because  the  coastal 
upwelling  phenomenon  is  usually  at  its  peak  during  that  month. 

B.  SPECIFIC  VOLUME  ANOMALY 
1.     Method 

Specific  volume  anomaly  (defined:  V  —  Vmp  where  V  is  specific  volume,  and 
subscripts  in  the  second  term  indicate  salinity  35  psu,  temperature  0°  C,  and  pressure  in 
db)  can  be  expressed  as  a  sum  of  constituents  related  to  the  parameters  that  define  the 
density  of  sea  water: 

S  =  (<5,  +  os  +  3st)  +  5[p  +  dsp  +  6stp. 

(1)     (2)     (3)       (4)      (5)      (6) 

Subscripts  t,  s,  and  p  stand  for  temperature,  salinity  and  pressure,  respectively.  In  this 
linear  construct,  specific  volume  anomaly,  also  termed  steric  anomaly,  is  the  sum  of  a 
part  due  to  variation  in  temperature  from  0°  C,  salinity  other  than  35  psu,  their  inter- 
action, and  the  interaction  of  deviant  temperature  and  salinity  with  pressure  (Sverdrup 


1   Unpublished  data  are  on  file  at  the  Department  of  Oceanography,  Naval  Postgraduate 
School,  Monterey,  CA. 


et  al.,  1942:57-59).  Since  the  term  5stp  is  much  smaller  than  the  other  terms,  it  is  virtually 
always  ignored,  and  is  not  considered  in  this  thesis.  The  sum  of  terms  in  parenthesis 
represents  the  total  combined  contribution  of  temperature  and  salinity  calculated  at 
surface  pressure,  is  called  the  thermosteric  anomaly,  Aff,  and  is  by  far  the  greatest  con- 
tribution to  the  total  specific  volume  anomaly.  Sverdrup  (1933)  introduced  the  use  of 
the  thermosteric  anomaly,  which  is  easily  determined  in  total  from  the  measured  or 
computed  surface  density,  as  a  computational  convenience.  After  the  density  anomaly 
calculated  at  zero  pressure  is  converted  to  A,„  only  the  two  values  d.p,  and  5!p  are  neces- 
sary to  determine  specific  volume  anomaly. 

Traditionally  then,  since  Sverdrup  (1933),  only  three  tables  (or  calculations)  are 
used  to  determine  d.  However,  it  can  sometimes  be  quite  useful  to  know  the  individual 
contributions  to  A„.  For  example,  if  <5S  is  small  or  constant  while  5e  is  large  and  fluctu- 
ating, a  dynamic  modeler  might  simplify  his  model  by  keeping  salinity  constant.  Simi- 
larly, if  a  feature  can  be  defined  by  its  salinity  signature,  an  observationalist  might 
identify  the  feature's  boundary  by  a  sharp  gradient  in  5S  . 

The  equation  of  state  for  seawater: 


'*        (1-Plks[p)' 


(4) 


where  k  is  the  secant  bulk  modulus  and  P  pressure,  translates  to  an  equation  for 
thermosteric  anomaly: 

Ast  =  -7 9.7266204  x  10~4  (5) 

51         PstO 


-  9.7266204  x  10    ,        (6) 


(Pw  =  ^O  +  s(£v")  +  s3/2£/,/) 


+  dnS* 


where  pw  is  the  density  of  reference  pure  water,  a,  b,  and  c  are  coefficients  in  power  se- 
ries on  temperature,  and  d0  is  a  coefficient  on  salinity  (UNESCO,  1987).  This  inverted 
sum  of  terms  in  S  and  T,  cannot  readily  be  analytically  sorted  into  parts  due  solely  to 
temperature,  salinity,  and  their  interaction.  Therefore,  some  other  means  of  defining 
6„  Ss,  and  5„  is  needed.    In  UNESCO  (1987)  equation  (5)  is  modeled  by: 


dstP  =  dsco  +  {dl5tp  -  d25l0)  +  {d,np  -  ds00) 

(I)         (2)         (3)  (4)        (5)  (7) 

(a)  (b)  (c) 

(Here  d  is  used  in  lieu  of  d  to  alleviate  confusion  that  might  otherwise  result  when 
comparing  terms  of  equation  (3)  to  terms  of  equation  (7)).  Clearly  specific  volume  cal- 
culated for  a  given  salinity  and  temperature  at  surface  pressure,  d:tQ,  is  precisely  the 
thermosteric  anomaly.  Term  (b),  the  difference  between  d  calculated  at  temperature  and 
depth  with  salinity  35  psu  and  d  calculated  for  the  same  temperature  but  at  surface 
pressure,  describes  the  interplay  of  temperature  and  pressure,  <),,,  .  Similarly,  term  (c) 
represents  Ssp.  (Terms  a,  b,  and  c  correspond  to  UNESCO,  1987,  tables  III,  V,  and  VI.) 
Now,  consistent  with  the  currently  accepted  equation  of  state  for  seawater 
(UNESCO,  1981),  how  does  one  model  the  individual  constituents  of  the  thermosteric 
anomaly?  Term  (3)  in  equation  (7)  is  the  value  of  S  that  would  result  at  zero  pressure 
if  salinity  were  a  constant  35  psu.  Therefore  [dstp  —  4sio)  represents  the  part  of  d  due  to 
salinity  deviation  from  35  psu  and  interaction  of  temperature  with  that  deviation  at  a 
certain  pressure  p  ,  or: 

dstp  ~  4?5r0  =  (^  +  «5,r  +  Stp  +  $sP)-  (8) 

Similarly: 

ds[p  ~  4oo  =  for  +  dtP  +  5sp  +  St).  (9) 

Subtraction  of  equation  (9)  from  equation  (S)  yields: 

^-Ss  =  d25[0-dm.  (10) 

It  would  appear  diSr0  and  ds00  can  serve  as  models  for  3,  and  Ss.  Using  that  assumption, 
equation  (3)  can  be  reconfigured  from  appropriate  elements  of  equation  (7): 


<5f  =  ^35/0 

$st  =  Aj,  ~  <>r  _  &s  =  </Jz0  -  ^35rt)  -  ^00 

*tp  =  d35tp-d35tO  (ll) 

^  =  dS0p  ~  "i00 
<^7/?  =  d35tO  +  ^00  +  (^rrO  ~~  ^35zO  ~~  4oo)  +  (^35tp  ~  ^35ri))  +  (^rO/>  —  ^jOo)1 

(1)       (2)  (3)  (4)  (5) 

Collection  of  terms  shows  equation  (1 1)  to  be  identical  to  equation  (7)  and  therefore  this 
model  to  be  consistent  with  the  UNESCO  (1987)  table  formulation.  Descriptively  this 
formulation  should  also  be  acceptable:  the  temperature  contribution  is  defined  as  the 
specific  volume  that  would  result  from  the  observed  temperature,  if  salinity  were  35  psu 
and  pressure  0  db;  the  salinity  contribution  at  the  same  depth  is  similarly  the  specific 
volume  anomaly  that  would  result  from  actual  salinity  and  0°  C;  the  combined  term  is 
simply  the  difference  of  the  sum  of  those  two  terms  and  Ayf  or  the  amount  by  which  their 
sum  misrepresents  AJf. 

It  follows  that  (0,  +  5tp)  can  be  represented  by  d35tp,  and  (<5,  +  Ssp)  by  ds<)p.  It  was 
decided  to  use  this  combined  temperature  and  temperature-pressure  term  in  order  to 
compare  total  saline  effect  to  temperature  effect  at  a  given  pressure.  As  a  notational 
convenience,  throughout  the  remainder  of  this  thesis  "dT"  and  "ds  "  will  be  used  in  lieu 
of  "(<),  +  c5J  "  and  "(Ss  +  5sp)  "  and  be  considered  synonymous  with  "dlUp "  and  "ds0p."  In 
this  formulation,  the  term  dsr  =  A„  -  ST-  Ss  will  include  pressure  effects.  Since  units 
used  to  describe  specific  volume  anomaly  will  in  every  case  be  \0~s  m3  kg~l,  they  will 
hereafter  be  deleted  in  the  text. 
2.     Results 

For  each  month  and  the  depths  10,  100,  200,  and  500  m,  each  of  the  fields  ST, 
6S,  and  ds[  were  computed  and  plotted  (see  note  1).  Generally,  the  specific  volume 
anomaly  field  consists  of  a  dome  centered  near  45°  N,  124°  W,  a  depression  centered 
near  San  Francisco,  and  a  ridge  parallel  and  adjacent  to  the  Baja  California  peninsula 
with  axis  south  to  north  about  3°  longitude  from  shore.  In  July  (Figure  6)  all  three 
features  are  well  established.  The  Columbia  River  outflow  is  at  its  peak  during  late  June 
and  early  July  (Huyer,  1983),  and  the  warm,  fresh  surface  plume  results  in  enhancement 
of  the  northern  dome.  July  is  also  a  period  of  strong  coastal  upwelling  (Huyer,  1983), 
and  the  cold,  salty  upwelled  water  results  in  intensification  of  the  central  depression. 
The  horizontal  gradient  in  d  decreases  with  depth,  particularly  in  the  northern  region. 


At  100  m  (Figure  7)  for  example,  there  is  no  evidence  of  the  warm,  fresh  core  so  pre- 
valent on  the  surface.  At  500  m  (Figure  8)  the  gradient  is  very  weak  everywhere  with 
dense  cores  centered  about  2°  off  the  coast  at  San  Francisco  and  the  northern  Baja 
California  peninsula.  An  important  observation  is  that  5st  (Figure  9)  values  are  one  to 
two  orders  of  magnitude  less  than  the  sum  of  dT  (Figure  10)  and  5S  (Figure  11).  This 
shows  the  steric  anomaly  to  be  very  nearly  equal  to  the  sum  of  a  temperature  and 
salinity  part  so  that  this  formulation  is  an  effective  linearization.  In  all  cases  <5;.  is  neg- 
ative, increasing  slightly  in  absolute  value  with  depth;  values  of -3  to  -6  being  common 
at  10  m  and  -6  to  -9  more  typical  at  500  m  (Figure  12).  Gradients  in  S!t  are  everywhere 
small;  however  the  cross  term  becomes  more  negative  in  areas  where  temperature  and 
pressure  act  together  to  increase  S.  and  less  negative  when  temperature  and  salinity  act 
in  conjunction  to  decrease  S.  When  salinity  and  temperature  oppose  change  in  specific 
volume  anomaly,  5st  remains  nearly  constant;  the  temperature-pressure  interaction  acts 
to  mitigate  change  in  specific  volume  anomaly  that  results  from  changes  in  either. 

Advection  of  subarctic  water  south  from  the  point  where  the  North  Pacific 
Current  feeds  the  CC,  results  in  cold  (low  Sr)  and  fresh  (high  ds)  tongue-like  structures 
superimposed  on  the  other  characteristic  features  of  the  CCS  described  in  part  below. 
In  the  immediate  vicinity  of  the  coast,  salinity  and  temperature  act  together  to  decrease 
steric  anomaly.  Near  the  surface  this  is  most  apparent  just  north  of  San  Francisco  Bay 
in  July  where  large  scale  upwelling  results  in  cold,  salty  water.  The  effects  on  density 
of  temperature  and  salinity  are  positively  correlated  with  specific  volume  anomaly  near 
45°  N,  127°  W,  where  the  warm,  fresh  discharge  of  the  Columbia  River  results  in  the 
geopotential  dome  described  above.  This  feature  is  strongest  in  summer,  weakest  in 
winter,  and  most  apparent  in  the  ds  field.  Such  seasonality  is  in  keeping  with  the 
Columbia  River's  outflow  cycle.  There  also  persists  a  fresh  dome  at  the  mouth  of  the 
Strait  of  Juan  de  Fuca,  but  it  is  restricted  to  the  immediate  vicinity  of  the  coast.  (Since 
the  fine  scale  salinity-temperature  interplay  associated  with  the  drainage  from 
Vancouver  Island  is  not  resolvable  on  the  coarse  scale  used  here,  activity  in  that  area  is 
generally  not  included  in  discussion  here.) 

Near  the  surface,  throughout  the  year,  mean  gradients  of  Sr  across  the  CCS 
(neglecting  small  scale  features  like  those  described  above)  are  at  least  50  percent  greater 
than  those  of  Ss.  At  a  depth  of  10  m,  absolute  values  of  Sr,  which  range  from  about  80 
to  400  north  to  south  through  the  year,  fluctuate  more  season  to  season  and  are  gener- 
ally larger  than  values  of  Ss  which  range  typically  from  50  to  250.  However  the  gradient 
in  <5r  diminishes    more  rapidly    with  depth    than  does  that  of  Ss  ;  the    result  is  that  at 
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Figure  6.      Steric  Anomaly,  July,  10  in:     Units  are   10  *  m3  kg  ',  and  the  contour 
interval  is  10  units. 
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Figure  7.      Steric  Anomaly,  July,  100  in:     Units  arc  lO-"  m3  kg\  and  the  contour 

interval  is  10  units. 
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Figure  8.      Steric  Anomaly,  July,  500  in:     Units  arc  10-8  m3  kg~l,  and  the  contour 
interval  is  10  units. 
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Figure  9.  Specific  Volume  Anomaly  (T-S),  July,  10  m:  Sst.  Units  arc 
1Q-8  mi  kg-1,  and  the  contour  interval  is  1  unit  (dashed  contour  interval 
is  0.5  units  from  solid  contours). 


Figure  10.      Specific  Volume  Anomaly  (T),  July,  10  m:     5r    Units  arc  10"  m5  Ag1, 
and  the  contour  interval  is  10  units. 
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I.      Specific  Volume  Anomaly  (S),  July,  10  m:     <55.    Units  are  10'*  m3  kg~l, 

and  the  contour  interval  is  10  units. 


Figure   12.        Specific    Volume    Anomaly    (T-S),    July,    500    in:     <5rt.       Units    are 

10-8  m3  kg-\  and  the  contour  interval  is  1  unit  (dashed  contour  interval 
is  0.5  units  from  solid  contours). 
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500  m.  3T  (Figure  13)  and  5S  (Figure  14)  fields  are  virtually  mirror  images  of  each  other. 
At  that  depth,  ranges  of  40  to  70  are  typical  of  5S,  as  compared  to  45  to  75  for  ST.  At 
an  intermediate  depth  of  100  m.  ranges  from  65  to  190  and  70  to  150  correspond  to  3T 
(Figure  15)  and  Ss  (Figure  16),  respectively,  while  ranges  of  60  to  135  and  40  to  85  are 
representative  of  5T  (Figure  17)  and  5S  (Figure  18)  at  200  m.  The  salinity  signature  of 
the  Columbia  River  is  practically  gone  at  100  m.  In  fact,  in  the  region  north  of  43°  N, 
from  100  m  down,  gradients  in  both  the  or  and  3S  fields  are  with  one  exception  consid- 
erably weaker  than  in  southern  areas.  The  exception  is  a  Ss  maximum  at  200  m  during 
October  in  the  data  rich  area  near  44°  N.  127°  W  (Figure  19).  This  unexplained,  but 
possibly  real  (the  central  value  is  based  on  22  observations  although  surrounding  mean 
values  are  from  as  few  as  10)  mean  feature  supports  a  gradient  of  8.2  x  10-'  m1  kg~[  in 
one  degree  of  longitude,  which  is  comparable  to  the  strongest  gradients  anywhere  in  the 
CCS  at  that  depth  and  much  stronger  than  any  other  in  the  northern  sector. 

At  100  m,  and  to  a  lesser  extent  200  m,  a  meridional  tongue  of  Sr  with  its  axis 
just  offshore  is  apparent.  That  is,  there  is  a  pronounced  offshore  cold  core  adjacent  to 
the  California  coast.  Cyclonic  flow  around  this  feature  is  consistent  with  nearshore 
subsurface  poleward  flow  typical  of  the  CCS. 

C.     VARIANCE  AND  CORRELATION  OF  THE  SPECIFIC  VOLUME  ANOMALY 
CONSTITUENTS 
1.     Method 

Often  the  absolute  value  of  a  quantity  is  of  less  importance  than  the  nature  of 
the  temporal  or  spatial  derivatives  of  that  quantity.  Statistics  generated  from  a  signif- 
icant number  of  observations  can  provide  information  concerning  those  differentials. 
Deviation  from  local  mean  values  can  cause  changes  in  local  gradient  or  prevent  change 
in  local  gradient  when  neighboring  values  change.  In  a  flow  supported  by  tight  gradi- 
ents, small  meanders  result  in  relatively  large  changes  in  local  values.  The  actual  vari- 
ance2  of  a  field  is  a  descriptor  of  likely  deviation;  therefore,  in  a  data  block,  large 
variance  may  indicate  large  gradient  or  large  variability  in  the  local  gradient.  The  precise 
effect  on  instantaneous  gradient,  of  course,  depends  on  behavior  of,  and  correlation 
with,  neighboring  data  blocks.  For  example,  if  all  nearby  mean  values  change  in  tandem 
with  each  other,  the  gradient  would  remain  constant  and  "no  change"  in  gradient  would 
result  from  variance.   Without  thorough  knowledge  of  the  nature  of  the  field  of  interest. 


2    Whether  or  not  the  sample  variances  calculated  here  are  valid  representations  of  true  vari- 
ance in  the  various  fields  is  a  relevant  question.   Section  IVF 1  addresses  this  point  in  detail. 
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Figure  13.  Specific  Volume  Anomaly  (T),  July,  500  m:  <5r.  Units  arc 
\Q-S  m3  kg~\  and  the  contour  interval  is  10  units  (dashed  contour  inter- 
val is  5  units  from  solid  contours). 
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Figure  14.  Specific  Volume  Anomaly  (S),  July,  500  in:  Ss.  Units  arc 
lO-8  m3  kg-1,  and  the  contour  interval  is  10  units  (dashed  contour  inter- 
val is  5  units  from  solid  contours). 
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Figure   15.        Specific    Volume    Anomaly    (T),    July,     100    in:     S,.        Units    arc 

10  R  /;;'  kg -'.  and  the  contour  interval  is  10  units. 


3V  ANOMALY  DUE  TO 
JAUNITY 
00  MPTflRS 
•:    U  M3/KG 


131)  W  120  W   126  W   124  W   122  W   120  W   11(1  W   llf-  W   111  W   112  W 

Figure   16.        Specific    Volume    Anomaly    (S),    July,     1 00    m:     Ss.        Units    arc 
I0"8  ms  kg~l,  and  the  contour  interval  is  10  units. 
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Figure   17.        Specific    Volume    Anomaly    (T),    July,    200    in:     <)7,       Units    arc 
10"  //;'  kg~l,  and  the  contour  interval  is  1()  units. 
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Figure    18.        Specific    Volume    Anomaly    (S),    July,    200    in:     Ss.        Units    arc 
I(H  m*  kg-1,  and  the  contour  interval  is  10  units. 
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it  is  impossible  to  quantify  gradient  from  sample  variance  calculation,  and  the  fact  that 
two  data  blocks  have  identical  variance  does  not  mean  that  they  have  identical  gradient. 
However,  determination  of  the  contribution  of  various  constituents  to  the  total  variance 
provides  strong  indication  of  the  contribution  of  those  constituents  to  the  development 
of  the  field  of  interest.  It  is  with  this  caveat  that  the  following  analysis  technique  is  used. 
Variance  of  the  specific  volume  anomaly  field  results  from  the  density  changes 
that  can  be  key  to  momentum  balances.  To  evaluate  the  importance  of  saline  variance 
to  specific  volume  anomaly  variance,  o}s  (variance  in  Ss)  should  be  compared  to  a] 
(variance  in  <3).  Because  equation  (3)  is  linear,  it  is  tempting  to  try  to  represent  the 
variance  of  6  as  the  sum  o^  the  variances  of  the  constituents.  If  this  were  possible, 
comparison  of  o]T  (variance  in  5T)  to  a],  for  example,  would  show  at  a  glance  how  much 
of  the  total  variance  is  due  to  temperature  variance  alone.  That  is,  whenever  the  ratio 
o\Jg]  remains  constant  and  large,  approaching  one  as  an  ideal  limit,  salinity  variance 
could  be  deemed  unimportant.  Conversely  if  a\T\a\  varies  significantly  over  time  or 
space  intervals  of  interest,  salinity  variance  could  be  considered  important.  However, 
because  the  constituents  of  specific  volume  anomaly  are  not  independent  of  each  other, 
this  approach  is  not  valid.  The  covariances  play  an  essential  role  in  determining  the 
variance  of  the  specific  volume  anomaly  field.  Therefore,  to  evaluate  a  term  on  the  right 
hand  side  of  equation  (3)  as  a  potential  estimator  of  S,  in  addition  to  variance  compar- 
ison, it  is  necessary  to  determine  how  that  term  correlates  with  3.  If.  for  example,  ST  is 
highly  correlated  with  6,  then  most  of  the  variance  in  <5  can  be  accounted  for  by  a  linear 
relationship  with  5T.  If  also  the  variance  of  Sr  is  of  similar  magnitude  to  that  of  5,  then 
ST  can  be  considered  a  good  estimator.  After  all,  if  6  changes  when  ST  changes,  in  similar 
amounts  and  directions,  then  ST  would  be  expected  to  model  S  well.  If  that  is  the  case, 
a  similar  look  at  the  behavior  of  6S  will  probably  show  it  contributes  little  to  total  vari- 
ance, and  salinity  might  be  considered  unimportant  to  the  development  of  the  specific 
volume  field. 

The  correlation  coefficient  is  useful  as  an  indicator  of  the  degree  of  correlation 
but  more  useful  as  an  indicator  of  the  sign  of  the  correlation.  The  coefficient  of  deter- 
mination, on  the  other  hand,  tells  exactly  how  much  of  the  variation  in  one  variable  is 
a  result  of  a  linear  relationship  with  the  other  (e.g.,  Walpole  and  Myers,  1985). 

These  ideas  can  be  used  to  evaluate  the  validity  of  the  use  of  temperature  alone 
to  determine  density.  Wherever  ST  is  highly  correlated  with  S.  and  has  variance  of  sim- 
ilar magnitude,  modeling  S  with  ST  (or  p  with  p{T))  should  give  acceptable  results.  If 
ST  is  highly  correlated  and  its  contribution  to  total  variance  nearly  constant,  even  if  not 


precisely  the  same  magnitude  as  total  variance,  a  ST  density  model  might  still  be  ac- 
ceptable because  it  will  predict  a  representative  gradient  that  will  differ  from  the  true 
gradient  by  a  nearly  constant  factor. 
2.     Results 

Variances  of  S,  5T,  and  5S  were  computed  for  each  level  of  each  data  block,  as 
were  the  covariances,  aSji  (covariance  between  ST  and  S)  and  ai$s  (covariance  between 
ds  and  <5  ).  Correlation  coefficients  between  S  and  Sr  (rss)  and  <5  and  os  (rss)  were  also 
computed: 


(12) 


where   x   represents  either  Ss  or  ST.    In  addition,  coefficients  of  determination,  r2ss,  and 
rlis  (as  used  in  Figures  22-29)  were  calculated. 

Maps  of  determination  coefficients  (given  the  sign  of  the  correlation  coefficient) 
for  all  months  and  the  depths  10,  100.  200  and  500  m  were  evaluated  (see  note  1). 
Multiplication  of  the  variance  in  ST  ,  by  r]r,  determines  how  much  of  the  total  variance 
in  5  is  due  to  variance  in  5T.  Below  about  100  m.  variance  in  the  5  field  is  very  small, 
too  small  to  allow  confidence  of  statistical  significance.  Analysis  of  ST  and  5S  corre- 
lations with  5  is  therefore  limited  to  the  upper  100  m.  Generally  5T  (Figure  20, 
Figure  21)  is  more  highly  correlated  with  5  than  is  5S.  Correlation  between  5S  and  5 
(Figure  22,  Figure  23)  is  usually  quite  weak  on  the  surface  and  stronger,  although  not 
as  large  as  r\bj  ,  at  100  m.  These  statements  are  true  only  in  a  very  general  sense.  For 
example,  frequently  in  both  fields  very  high  or  low  correlations  exist  in  the  immediate 
vicinity  of  each  other.  Another  exception  is  that,  as  expected,  immediately  offshore  of 
a  freshwater  source,  dT  correlations  are  occasionally  weak,  or  nonexistent.  Usually  when 
the  <5r  correlation  is  weak,  the  5S  correlation  is  very  high.  For  example  at  10  m  in  Feb- 
ruary at  the  mouth  of  San  Francisco  Bay,  r]iT  (Figure  24)  is  0.0  while  r]Ss  (Figure  25)  is 
1.0.  Both  correlations  can  be  very  high,  predictably,  for  instance,  in  the  Columbia  River 
plume,  where  warm,  fresh  water  results  in  high  specific  volume  anomaly.  On  the  other 
hand,  the  opposite  situation  exists  well  offshore  of  Baja  California,  near  Guadalupe  Is- 
land. There,  in  late  winter  and  spring,  very  low  6T  correlations  (Figure  26)  coexist  with 
very  low  r\iS  values  (Figure  27).    Variance  in  both  6S  and  Sr  is  very  high  in  relation  to 


5,  but  the  two  are  so  strongly  in  opposition  that  neither  correlates  well  with  the  total 
variance. 

Throughout  the  year,  within  a  two  or  three  degree  radius  of  31°  X,  119°  W  in 
the  upper  100  m,  r]6j.  (Figure  28)  is  consistently  near  0.9  and  the  variance  in  5T 
(Figure  29)  accounts  for  a  reasonably  consistent  part  of  the  variance  in  <5.  Perhaps  in 
this  area,  the  assumption  of  constant  salinity  would  introduce  minimal  error  in  a  density 
model. 

D.     MEAN  DIFFERENCES  OF  THE  SPECIFIC  VOLUME  ANOMALY 
CONSTITUENTS 

1.     Method 

If  an  observed  quantity,  x.  is  modeled  by  a  predictor  x\  then  (x  —  xK)  can  be 
considered  the  error  associated  with  a  single  observation.  If  that  error  is  squared  to 
eliminate  sign  bias  and  the  sum  of  such  squared  error  averaged,  the  square  root  of  the 
result  can  be  used  as  an  indicator  of  the  expected  deviation  and  is  termed  the  Root  Mean 
Square  Error  (RMSE).  The  R.V1SE  normalized  by  the  mean  value  of  x  then  provides 
an  estimate  of  the  magnitude  of  the  expected  error  in  terms  of  percent  of  the  "correct" 
value.    If  3  is  the  observed  quantity  and  Sx  the  predictor: 


Z 


-  dxh 


I") 

=RMSE„,  (13) 


If  for  example  x  is  T,  the  result,  RMSE(T),  provides  an  indication  of  how  well  5T  ap- 
proximates 6  in  a  particular  data  block. 
2.     Results 

Normalized  RMSE  values  were  calculated  for  5r  and  Ss.  As  discussed  in  section 
IVB,  both  dT  and  5S  are  everywhere  greater  than  zero  and  each  greater  than  the  absolute 
value  of  the  sum  of  the  remaining  specific  volume  anomaly  constituents.  Therefore  5T 
and  5S  are  necessarily  less  than  6,  and  the  sign  of  the  expected  error  is  known  (the  "root" 
and  "squared"  point  of  the  RMSE  technique  is  moot,  and  RMSE  is  equivalent  to  mean 
difference).   The  specific  volume  anomaly  will  always  be  underestimated. 

As  discussed  earlier,  it  is  the  gradient  in  the  error  field  that  is  likely  more  sig- 
nificant than  magnitude  of  a  single  error.  After  all,  use  of  the  actual  mean  salinity  rather 
than  standard  35  psu  can  remove  a  constant  offset.    The  plots  of  RMSE(T)  provide  a 
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Figure  20.      Specific  Volume  Coefficient  of  Determination  (  I),  July,  10  in:     r\if 
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Figure  21.      Specific  Volume  Coefficient  of  Determination  (T),  July,  100  in:     r\Sr 
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Figure  22.      Specific  Volume  Coefficient  of  Determination  (S),  July,  10  in: 
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Figure  23.      Specific  Volume  Coefficient  of  Determination  (S),  July,  100  in: 
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Figure  24.      Specific  Volume  Coefficient  of  Determination  (T),  February,  10  in: 
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Figure  25.       Specific  Volume  Coefficient  of  Determination  (S),  February,  10  m: 
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Figure  26.      Specific  Volume  Coefficient  of  Determination  ( T),  March,  10  m:     r}Sr. 
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Figure  27.      Specific  Volume  Coefficient  of  Determination  (S),  March,  10  in: 
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Figure  28.      Specific  Volume  Coefficient  of  Determination  (T),  January,  10  m:     r\if 
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Figure  29.        Variance   in  Steric   Anomaly   from   Temperature,   January,    10    in: 
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clear  large  scale  picture  of  the  price  of  neglecting  salinity  variation  in  specific  volume 
anomaly  calculations  (e.g..  Figure  30).  Near  the  surface  that  neglect  would  result  in  a 
40  to  50  percent  underestimation  in  the  north-south  gradient,  and  a  10  percent  under- 
estimate west  to  east,  with  the  largest  errors  tending  to  occur  in  winter  months.  The 
meridional  gradient  persists  but  decreases  to  a  year-round  20  percent  at  500  m 
(Figure  31). 

Because  ST  and  Ss  are  by  far  the  greatest  contributors  to  total  6,  data  block 
values  of  RMSE(S)  added  to  the  corresponding  values  for  dT  usually  are  nearly  equal  to 
one  (compare  Figure  32  to  Figure  30).  Nevertheless,  Ss  mean  difference  plots  are  il- 
lustrative of  the  part  5S  plays  in  determining  total  steric  anomaly.  The  large  scale  op- 
position of  gradients  in  ds  and  ST,  in  virtually  every  direction,  is  clear  from  analysis  of 
these  distributions.  At  500  m,  in  the  southern  section  of  the  CCS  there  is  a  pronounced 
increase  in  the  zonal  gradient  in  the  Ss  RMSE  field  (Figure  33).  salinity  making  a  larger 
contribution  to  <5  offshore  than  nearshore.  Some  smaller  scale  complexities  are  also 
discernible.  For  instance,  near  the  surface,  across  the  mouth  of  San  Francisco  Bay  in 
September  (Figure  34),  there  is  a  twenty  percent  change  in  the  ds  contribution  to  3,  ac- 
companied by  a  very  flat  gradient  in  RMSE  for  dT.  This  phenomenon  is  illustrative  of 
the  local  nature  of  saline  variability,  and  temperature-salinity  interaction  in  the  CCS. 

An  interesting  result  of  this  analysis  is  that  there  are  large  areas,  virtually  any 
six  degree  latitude  by  four  degree  longitude  block,  where  the  assumption  of  no  salinity 
variation  would  apparently  result  in  less  than  a  ten  percent  error  in  the  resultant  specific 
volume  gradient.  That  is  particularly  true  if  areas  immediately  adjacent  to  the  coast  are 
excluded. 

E.      DYNAMIC   HEIGHT 

1.     Method 

Dynamic  height,  O,  is  a  measure  of  how  the  density  of  a  column  of  water  differs 
from  the  density  of  a  standard  column  of  the  same  depth  and  is  calculated  by  integrating 
the  specific  volume  anomaly  over  depth  from  an  assumed  layer  of  no  motion  to  a  depth 
of  interest.  It  is  the  horizontal  gradient  in  dynamic  height  that  determines  the 
geostrophic  velocity  field.  If  the  contributions  to  that  gradient  from  salinity  are  signif- 
icant, they  should  be  considered  when  momentum  calculations  are  performed.  If  the 
mean  differences  (d  -  ST),  where  (  )  denotes  the  mean,  in  columns  are  integrated  with 
respect  to  pressure,  the  resulting  horizontal  distribution  is  a  good  measure  of  the  effect 
of  salinity  on  the  geostrophic  velocity  field.    That  is,  contours  of  cj>  -  Oiirp  =  O1  define 
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Figure  30.      Specific  Volume  Anomaly  RMSE  (T),  July,  10  in. 
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Figure  31.      Specific  Volume  Anomaly  RMSE  (T),  July,  500  in. 
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Figure  32.      Specific  Volume  Anomaly  RMSE  (S),  July,  10  in. 
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Figure  33.      Specific  Volume  Anomaly  RIY1SE  (S),  July,  500  in. 
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Figure  3-4.      Specific  Volume  Anomaly  RMSE  (S),  September,  10  m. 


mean  geostrophic  flow  that  will  not  be  accounted  for  if  5T  is  used  as  an  estimator  of  <5. 
Therefore,  if  gradients  in  the  difference  field  are  large,  salinity  is  important. 
2.     Results 

Plots  of  monthly  mean  values  of®,  ®35f?,  and  <D\  obtained  by  integrating  from 
an  assumed  layer  of  no  motion  at  500  m,  to  200  m  (200,500),  100  m  (100:500),  and  the 
surface  (0/500)  were  plotted  and  evaluated.  The  d>  results,  where  they  overlap,  are  in 
very  close  agreement  with  those  of  Wyllie  (1966).  Standard  deviations  are  in  virtually 
every  case  considerably  more  than  an  order  of  magnitude  less  than  mean  values.  Even 
so,  this  small  variance  still  indicates  that  contour  locations  are  doubtful  by  more  than 
the  width  of  a  contour  interval  on  these  maps.  Therefore,  any  conclusions  drawn  con- 
cerning small  scale  features  are  suspect.  It  should  also  be  remembered  that  these  are 
maps  of  mean  fields,  representing  climatology  rather  than  instantaneous  flows. 

Examination  of  the  0  500  dynamic  height  field  of  the  CCS  (e.g.,  Figure  35  and 
Figure  36).  leads  to  division  of  the  area  into  three  regions.  North  of  43°  N  the  diver- 
gence of  the  North  Pacific  Current  and  freshwater  discharges  from  the  Columbia  River 
and  Strait  of  Juan  de  Fuca  result  in  disorganized  mean  geostrophic  gradients.  A  split 
of  the  flow  toward  the  north  and  south  is  superimposed  on  weak  circulations  around 
warm,  freshwater  plumes.  The  area  along  the  U.S.  coast  from  30°  to  3S°  N  is  dominated 
by  cyclonic  geostrophic  circulation  inshore,  and  flow  generally  to  the  southeast  offshore. 
The  strength  of  the  southward  flow  and  degree  of  cyclonic  turning  onshore  is  seasonal. 
Southward  flow  is  strongest  in  summer  and  weakest  in  winter.  Conversely,  onshore 
cyclonic  turning  is  least  in  summer,  greatest  in  winter.  There  is  another  cyclonic  center 
just  oilshore  of  Baja  California  near  25°  N.  Again  the  cyclonic  departure  from  a  gen- 
erally southeastward  flow  is  seasonal  as  described  above.  For  each  month  the  dynamic 
height  difference  fields  (0\)  (Figure  37)  show  inclusion  of  horizontal  salinity  gradients 
to  be  responsible  for  a  nearshore,  equatorward  geostrophic  component  from  San 
Francisco  to  San  Diego  and  a  southwestward  flow  along  the  coast  of  Baja  California. 
This  predictable  geostrophic  pattern  is  due  to  the  persistent  salinity  pattern:  fresh  to 
saline,  offshore  toward  shore,  with  a  more  meridional,  fresh  to  salty  north  to  south, 
gradient  in  the  southern  section.  What  is  even  more  interesting  is  the  fact  that 
throughout  the  year,  the  magnitude  of  the  anomalous  flow  implied  by  the  <1>V  field  is  on 
the  same  order  as  flow  associated  with  V®  itself,  and  perhaps  greater  in  some  areas  than 
the  magnitude  of  the  flow  predicted  by  V<F>35r;.  (compare  figures  36-3S). 

In  the  northern  region,  in  each  month.  V®35r„  fails  to  account  for  a  persistent 
south-southwestward  component  from  the  mouth  of  the  Columbia  River  that  varies  in 
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Figure  35.      Dynamic  Height,  0/500  m,  January:     <I\  the  contour  interval  is  0.2  and 
units  arc   dynamic  decimeters. 
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Figure  36.       Dynamic  Height,  0/500  m,  July:     <I\  the  contour  interval  is  0.2  and 
units  are   dvnamic  decimeters. 
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Figure  37.       Dynamic  Height  Difference,  0/500  in,  July:     <I>%  =  <I>  -  <!>„„,  the  con- 
tour interval  is  0.2  and  units  are   dynamic  decimeters. 
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Figure  3<S.      Dynamic  Height  (Salinity  35),  0/500  in,  July:     <!>,<„.  the  contour  inter- 
val is  0.2  and  units  arc   dynamic  decimeters. 


intensity  with  the  magnitude  of  discharge  from  that  river.  However,  seaward  of  that 
southward  component  the  difference  field  is  very  slack.  That  is.  the  dynamic  height  field 
estimated  from  temperature  gradient  alone  very  nearly  accounts  for  the  geostrophic  flow 
pattern  in  the  region  where  surface  salinity  has  the  largest  annual  variance.  This  ap- 
parent contradiction  is  explained  by  the  fact  that  the  horizontal  salinity  gradient  here 
decreases  rapidly  with  depth.  At  depth,  throughout  the  year,  there  is  virtually  no  sig- 
nificant salinity  gradient  in  this  northern  area.  At  200  m  for  example,  over  a  distance 
where  mean  temperature  may  van-  by  0.4°  C,  salinity  changes  less  than  0.1  psu.  In  the 
southern  area,  by  contrast,  a  0.5°  C  temperature  change  is  accompanied  typically  by  a 
change  of  0.3  psu.  The  result  is  that  in  the  north,  the  depth  integrated  effect  of  salinity 
on  geostrophic  flow  only  approaches  that  of  temperature  in  the  immediate  vicinity  of  the 
sharpest  upper  level  salinity  gradient. 

Analyses  of  dynamic  height  and  geostrophic  flow  at  other  levels  (e.g.,  200500 
and  100  500  m.  Figures  39-44)  show  relative  gradients  very  similar  to  those  previously 
discussed.  VOv  is  of  similar  order  as  VO  and  V<friSlp  at  all  depths,  and  the  general  pattern 
of  the  difference  field  persists  with  depth.  North  of  northern  Baja  California,  the  differ- 
ence field  shows  a  strong  offshore-onshore  gradient  indicative  of  higher  salinities  inshore 
in  the  poleward  flow.  This  result  is  consistent  with  the  idea  that  the  poleward  flow  can 
be  traced  by  its  high  salinity  (e.g.,  Lynn  and  Simpson.  1987).  That  is.  the  inshore 
poleward  undercurrent  may  carry  the  saline  signature  of  Equatorial  water  which  is  to  the 
south  of  the  CCS.  That  salinity  increases  significantly  toward  the  south  is  evident  from 
the  fact  that  south  of  about  34°  N  there  is  a  strong  alongshore  north  to  south  gradient 
in  the  difference  field.  The  southward  salinity  increase  is  also  apparent  in  the  salinity 
fields  (Figure  4).  The  100  500  and  200/500  m  dynamic  height  calculations  show  in- 
creased cyclonic  turning  associated  with  the  large  scale  California  coastal  features  de- 
scribed earlier.  That  is.  the  dense  cores  associated  with  what  look  like  "c"  shaped  troughs 
centered  at  the  California  coast  and  Baja  California  peninsula  on  the  0/500  m  plots  move 
offshore  with  depth  so  that  at  200  m,  and  possibly  even  at  100  m,  a  closed  geostrophic 
circulation  with  coastal  poleward  flow  exists.  Because  of  the  neglected  zonal  salinity 
gradient.  <t>Xv  greatly  overestimates  the  strength  and  horizontal  extent  of  the  poleward 
flow. 

Due  to  the  absence  of  the  freshwater  runoff  below  the  surface,  the  divergence 
of  the  North  Pacific  Current  is  relatively  clear  on  the  100'500  and  200  500  m  plots.  This 
divergence  area  seems  to  be  farther  south  in  the  winter  than  summer,  as  noted  by 
Pickard  and  Emery  (19S2).  Although  the  divergence  occurs  in  a  region  of  generally  poor 


data  coverage,  the  seasonal  migration  seems  to  be  from  near  40°  N  in  the  winter  to 
45°  N  in  summer,  rather  than  the  45°  N  to  50°  N  range  reported  by  Pickard  and  Emery 
(1982). 

The  focus  of  this  thesis  is  large  scale;  however,  using  the  rule  of  thumb  that 
features  whose  wavelengths  exceed  2  grid  points  in  extent  may  be  resolved,  it  may  be 
possible  to  use  these  calculations  to  gain  insight  into  some  relatively  small  scale  features. 
In  Appendix  B.  discussion  of  some  smaller  scale  effects  of  salinity  on  dynamic  height  in 
the  CCS  is  provided. 

F.     PROPAGATION  OF  T-S  VARIANCE  THROUGH  THE  EOS 

1.     Method 

An  alternate  approach  to  examine  the  relative  importance  of  constituent  vari- 
ance to  density  variance  is  to  propagate  the  variance  associated  with  each  of  the  inde- 
pendent variables  through  the  complete  nonlinear  equation  of  state.  By  the  chain  rule, 
if  density  is  strictly  a  function  of  salinity,  temperature,  and  pressure  then: 

dp  cp  en 

Ap  =  -4-At  +  -4-As  +  -r-Ap.  (14) 

CI  cs  cp 

where  (A)  indicates  departure  from  mean  values  and  p  is  mean  density.  Because  all  ob- 
servations are  on  standard  depth  levels  and  pressure  changes  on  those  levels  are  small, 
even  in  comparison  to  error  associated  with  depth  determination.  Ap  can  be  neglected. 
Squaring  both  sides  of  (14).  summing  over  a  number  of  observations,  and  dividing  by 
the  number  of  observations  less  one  (e.g.,  Greenwalt  and  Shultz,  1962)  leads  to  : 

2      t  ep  ^2  2  ,  ,  OP  v2  2^~,  dP  v  ?P  v  „« 

where  a2  represents  variance  (either  sample  or  population),  a  denotes  covariance.  and 
subscripts  represent  variables  to  which  the  statistic  applies.  As  demonstrated  by 
Lillibridge  (19S9),  the  equation  of  state,  EOS(80),  a  polynomial,  is  relatively  easy  to 
differentiate.  Therefore  the  contribution  to  c\  from  temperature,  salinity,  and  their 
interaction  is  easily  determined  over  any  interval  of  interest.  This  technique  can  be  as 
readily  applied  to  specific  volume  or  any  of  the  specific  volume  anomaly  constituents 
and  is  desirable  because  it  employs  the  full  nonlinear  equation  of  state.  Insofar  as  the 
statistics  used  are  accurate  estimators,  the  equation  of  state  itself  accurate,  and  the 
pressure  levels  constant,  the  expression  is  valid.    The  method  is  particularly  useful  when 
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Figure  39.      Dynamic  Height,  l()()/5()()  in,  July:     <I\  the  contour  interval  is  0.2  and 
units  are   dynamic  decimeters. 
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Figure  40.      Dynamic  Height  (Salinity  35),  100/500  in,  July:     O^,  the  contour  in- 
terval is  0.2  and  units  arc   dynamic  decimeters. 
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Figure  41.        Dynamic   Height    Difference,    100/500   m,   July:     d)'  =  <I>  -  Q>Xlp,   the 

contour  interval  is  0.2  and  units  arc   dvnamic  decimeters. 


DYNAMIC  HEIGHT 

200/500 

IN  DYNAMIC  DECIMETERS 

JULY 


133     TVHI     W  129     W  127     W125     W123     W121     W119     W117     W115     W  II 3     Will     W109     W  1  (V 


Figure  42.       Dynamic  Height,  200/500  in,  July:     <l\  the  contour  interval  is  0.2  and 
units  arc    dynamic  decimeters. 
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Figure 


43.      Dynamic  Height  (Salinity  35),  200/500  in,  July:     <I>35„,  the  contour  in- 
terval is  0.2  and  units  arc   dynamic  decimeters. 
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Figure  44.       Dynamic   Height    Difference,   200/500  m,   July:     $'  =  O  -  Q>„„  the 

contour  interval  is  0.2  and  units  arc   dynamic  decimeters. 


temperature  and  salinity  ranges  over  a  time-space  interval  of  interest  are  reasonably  well 
known  and  a  quick  assessment  of  the  relative  importance  of  constituent  variances  de- 
sired.   It  also  is  computationally  quite  efficient  because  the  equation  of  state  need  only 
be  invoked  once  to  determine  total  and  constituent  variance  over  a  sample  of  interest. 
Let  the  terms  on  the  right  hand  side  of  equation  (15)  be  represented  as: 


2         .,  cp       dp 

As  mentioned  earlier,  to  increase  the  number  of  observations  and  hence  statistical  va- 
lidity, these  data  block  variance  estimates  were  computed  over  three-month  seasons, 
with  winter  defined  as  January.  February.  March. 

Potential  energy  (PE)  of  a  column  of  water  can  be  determined  by  integrating 
PE  =  pgz.  where  g  is  the  acceleration  due  to  gravity,  considered  constant,  and  z  is  height 
above  some  reference  level,  over  depth  (e.g..  rofonofT,  1962).  It  follows  that  on  a  con- 
stant depth  level  (or  very  thin  layer)  p  is  the  only  variable  in  PE.  and  the  horizontal 
density  field  can  be  used  to  represent  potential  energy.  Change  in  the  density  field,  either 
over  time  or  distance,  then  represents  change  in  the  potential  energy  field.  Equivalently, 
variance  from  the  mean  density  field  can  be  considered  as  an  indication  of  energetic  ac- 
tivity or  potential.  In  general,  relatively  high  density  variance  should  correspond  to 
relatively  high  energetic  activity.  Therefore,  if  sample  variance  computed  over  some  time 
or  space  interval  is  an  accurate  estimator  of  the  true  density  variance,  sample  variance 
may  provide  indication  of  energetic  activities.  Dantzler  (1977)  used  similar  statistical 
techniques  applied  to  vertical  excursions  of  the  thermocline  to  calculate  North  Atlantic 
potential  energy  levels. 
2.     Results 

Maps  of  g]  and  percent  contributions  of  salinity.  {a2pJa]  x  100),  temperature, 
(alJo;x  100).  and  covariance.  (op  Ja2px  KM)),  were  evaluated.  Near  the  surface,  large 
density  variance  occurs  in  areas  of  maximum  fluctuation  in  precipitation,  heating,  or 
freshwater  influx  for  example.  Accordingly,  Figure  45  indicates  high  variance  near  the 
Columbia  River  outflow,  and  in  the  upwelling  region  along  the  coast  of  California.  As 
expected,  salinity  variance  is  dominant  in  the  Columbia  River  plume  (Figure  46).    Also 


as  anticipated,  in  the  upwelling  region,  temperature  and  salinity  effects  (compare 
Figure  47  and  Figure  46)  are  positively  correlated  and  of  similar  relative  importance  to 
total  density  variance.  These  near-surface  observations  are  consistent  with  results  al- 
ready presented  and  so  provide  evidence  that  this  technique  may  be  valid. 

Away  from  the  surface,  areas  of  highest  indicated  variance  may  be  dynamically 
active  regions  in  the  CCS.  Regions  of  maximum  a;  might  correspond  to  the  mean  po- 
sitions of  currents  or  eddies.  If  so.  equation  (15)  can  be  used  to  determine  the  relative 
contributions  of  temperature,  salinity,  and  their  interaction,  to  the  CCS  itself.  In  sup- 
port of  this  idea,  the  locations  of  maximum  a\  seem  to  more  closely  correspond  to  CC 
and  CU  locations  as  depth  increases.  For  example  areas  with  o\  at  least  twice  as  high 
as  neighboring  data  blocks  at  200  m  (Figure  48)  correspond  reasonably  well  with  the 
areas  of  tight  gradient  in  O  (200  500  m)  apparent  in  Figure  42.  In  the  regions  of  high 
variance,  salinity,  with  variance  ratios  [g-p)\g\  typically  from  50  to  200,  seems  to  con- 
tribute more  to  total  variance  than  temperature,  with  variance  ratios  {o2pJo2p)  in  the  high 
total  variance  region  of  5  to  150.  At  500  m,  salinity  contributes  virtually  all  of  the  density 
variance  (Figure  49)  in  a  region  that  resembles  the  characteristic  cyclonic  circulation  of 
the  CCS.  The  high  correlation  of  saline  contribution  to  density  variance  indicates  that 
the  density  gradient  is  largely  due  to  salinity  gradients  at  this  depth.  Figure  50  shows 
some  areas  at  500  m  with  total  density  variance  values  on  the  same  order  as  at  200  m 
and  these  areas  have  values  at  least  twice  as  high  as  neighboring  data  blocks.  The  lo- 
cation of  these  regions  match  the  location  of  some  of  the  energetic  portions  of  the  CCS 
discussed  above.  These  observations  provide  evidence  that  the  density  gradient  at 
500  m  may  be  significant  and  is  largely  induced  by  salinity  gradients.  This  speculation 
is  consistent  with  that  of  Lynn  and  Simpson  (1989). 

If  quarterly  mean  values  of  salinity  and  temperature  are  substituted  for  /  and  5 
in  equation  (15).  and  statistics  are  generated  from  these  mean  values  by  summing  across 
the  length  or  width  of  the  CCS,  equation  (16)  will  represent  contributions  of  salinity, 
temperature,  and  their  covariance  to  the  total  density  variance  and  may  indicate  their 
importance  to  the  mean  density  gradient.  Such  calculations  were  performed  to  obtain 
statistics  cross-shore  at  each  latitude  that  had  at  least  five  seasonal  means  of  at  least  10 
observations,  and  alongshore,  across  data  blocks  equidistant  from  the  coast. 
Alongshore  lines  are  numbered  1  to  8  from  inshore  out,  and  cross-shore  lines  are  iden- 
tified by  latitude  (Figure  51).  The  results  at  the  surface  and  100  m  were  plotted  as 
histograms.  Inshore,  alongshore  lines  1  (Figure  52)  and  2  (Figure  53)  show  large  year- 
round  ti  tal  variance  a:  the  surface.     In  the  northern  sector  where  line  1  crosses  the 
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Figure  45.      Total  Density  Variance,  10  in,  Summer:       Units  arc  (kg m  ')?  x  10\ 
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Figure  46.       Salinity  Contribution  to  Total  Density  Variance,   10  m,  Summer: 


VARIANCE  RATIO  {%) 
TEMPERATURE/rOTAL 
10  METERS 
JUL/AUG/SEP 


130    W   120    W       126    W       12-1    W       122    W       120    W       1 10    W       11G    W        1H     W        112    W 

Figure  47.         Temperature     Contribution     to     Total     Density     Variance,     10     in, 
Summer:      o\  \o\  x  100. 
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Figure  48.      Total  Density  Variance,  200  in,  Slimmer:      Units  arc  {kgm  3)J  x  I<)\ 
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Figure  49.       Salinity  Contribution  to  Total  Density  Variance,  500  m,  Summer: 
alja2px  100. 
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F  igure  50.      Total  Density  Variance,  500  in,  Summer:      Units  arc  (kg  nr3)2  x  I0\ 


Columbia  River  outflow  region,  the  greatest  surface  contribution  is  clearly  from  salinity. 
On  this  line  even  at  100  m,  salinity  variance  dominates  alongshore  except  in  winter  when 
temperature  variance  is  apparently  slightly  more  important.  Line  2,  which  is  influenced 
by  runoff  the  length  of  the  California  coast,  as  well  as  by  Columbia  River  outflow,  shows 
surface  salinity  as  most  important  in  fall,  and  a  strong  saline  signature  year  round.  At 
100  m,  temperature  is  clearly  dominant  in  the  low  runoff,  winter  season,  and  at  least  of 
equal  importance  the  rest  of  the  year.  Northern  cross-shore  sections  (only  sections  at 
46°.  47°,  49°  N  in  the  autumn,  and  26°-35°  N  met  the  minimum  number  of  observations 
criterion)  show  salinity  variance  as  key  during  the  fall  (e.g.,  Figure  54). 

Alongshore,  the  covariance  is  always  negative,  as  must  be  expected  when  the 
large  scale  opposing  meridional  trends  in  temperature  and  salinity  are  considered. 
Cross-shore,  except  in  the  northern  area  where  warm,  fresh  Columbia  River  outflow  re- 
sults in  strong  positive  interactions,  covariance  tends  to  be  less  important  and  of  either 
sign,  reflecting  the  variability  of  the  nearshore  waters  (e.g.,  Figure  55). 

Cross-shore  variance  at  100  m  is  very  small  in  the  north,  generally  an  order  of 
magnitude  less  than  the  0.1  to  0.5  {kgrrr*)2  typical  at  the  surface.  In  the  south,  cross- 
shore  variance  at  100  m  is  of  comparable  magnitude  to  that  at  the  surface. 

Histograms  generated  for  areas  where  relative  contributions  to  density  variance 
of  temperature,  salinity,  and  their  interaction  are  the  same  will  have  similarly  shaped 
histograms.  Therefore,  it  is  reasonable  to  expect  that,  in  the  same  water  mass,  areas 
undergoing  similar  energetic  processes  will  have  similarly  shaped  histograms.  This 
reasoning  implies  that  if  the  shape  of  the  histogram  for  a  particular  section  is  the  same 
at  the  surface  and  100  m.  the  two  levels  may  well  be  undergoing  the  same  dynamics 
and  or  thermodynamics;  that  is.  they  are  dynamically  coupled.  In  the  autumn,  the 
graphs  for  areas  south  of  31°  N  at  both  depths  are  very  similar  in  shape  (e.g., 
Figure  56),  indicating  such  coupling  may  exist.  Perhaps  surface  cooling  which  occurs 
during  the  months  of  October  through  December  results  in  deeper  mixed  layers  south 
of  31°  X. 

In  this  southern  region,  except  when  total  variance  is  very  low.  temperature 
provides  the  largest  contribution  to  cross-shore  gradients.  There  is  a  tendency  for 
greatest  temperature  variance,  and  therefore  total  variance,  to  occur  in  fall,  further  evi- 
dence that  surface  cooling  is  greatest  during  that  period. 

Alongshore  lines  also  show  a  noticeable,  if  less  pronounced,  increase  in  variance 
during  the  autumn.  During  that  season,  the  alongshore  sections  seaward  of  line  2  show 
the  surface  and  100  m  layers  to  be  well  coupled  the  length  of  the  CCS  (e.g..  Figure  57). 


Figure  51.  Alongshore  and  Cross-Shore  Density  Variance  Contribution:  Con- 
tributions to  total  density  variance  by  temperature,  salinity,  and  their 
interaction  were  calculated  alone  the  lines  shown. 
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Figure  52.  Alongshore  Line  1,  Contributions  to  Density  Variance:  Contribution 
to  total  density  variance  (labeled  T)  by  temperature  (t),  salinity  (s),  and 
their  interaction  (t/s)  arc  depicted  here  by  season  at  the  surface  and  100 
m.  Units  of  density  variance  here  arc  (kg  m~3)2. 


Figure  53.  Alongshore  Line  2,  Contributions  to  Density  Variance:  Contribution 
to  total  density  variance  (labeled  T)  by  temperature  (t),  salinity  (s),  and 
their  interaction  (t/s)  are  depicted  here  by  season  at  the  surface  and  100 
in.  Units  of  density  variance  here  arc  (kg  m  i)2. 
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Figure  54.  Cross-shore  Line  46,  Contributions  to  Density  Variance:  Contribution 
to  total  density  variance  (labeled  T)  by  temperature  (t),  salinity  (s),  and 
their  interaction  (t/s)  arc  depicted  here  by  season  at  the  surface  and  100 
m.  Units  of  density  variance  here  are  (kg  trr*)7. 


Figure  55.  Cross-shore  Line  26,  Contributions  to  Density  Variance:  Contribution 
to  total  density  variance  (labeled  T)  by  temperature  (t),  salinity  (s),  and 
their  interaction  (t/s)  arc  depicted  here  by  season  at  the  surface  and  100 
in.  Units  of  density  variance  here  are  (kg  m  -)J. 


Again  a  possible  interpretation  is  that  surface  cooling  in  the  fall  results  in  deeper  mixing. 
If  this  interpretation  is  valid  then  mixing  occurs  offshore  for  the  length  of  the  CCS. 

It  may  be  argued  that  over  any  areas  or  periods  where  the  relative  contributions 
to  total  variance  remain  constant,  a  density  model  based  on  temperature  alone  can  ac- 
curately represent  gradients,  erring  only  by  a  constant  factor.  In  that  case,  intersections 
of  a  group  of  adjacent  alongshore  sections  sharing  a  characteristic  histogram  shape,  with 
a  group  of  cross-shore  sections  that  have  a  common  histogram  shape,  will  highlight  such 
regions.  Subject  to  individual  interpretation,  some  seasonal  combinations  of  lines  2-5 
and  26-31  may  meet  the  criteria.  For  example,  the  autumn  histograms  for  lines  3 
(Figure  58)  and  4  (Figure  59)  are  nearly  the  same  shape,  as  are  the  autumn  graphs  of 
lines  27.  28,  and  29  (Figure  60,  Figure  56  and  Figure  61,  respectively).  The  implication 
here  is  that  during  October  through  December,  at  least  in  the  upper  100  m,  within  an 
area  2°-4c  from  shore,  between  27°  and  30°N,  a  temperature-only  density  model  would 
overestimate  the  alongshore  gradient  by  a  constant  small  factor  and  underestimate  the 
cross-shore  gradient  by  a  similar  amount.  Again  the  indication  is  that  offshore,  in  select 
three  or  four  degree  blocks,  salinity  variance  may  not  be  essential  to  an  accurate  model 
of  the  dynamic  field.  (Using  this  technique  to  generate  shorter  alongshore  cuts,  of 
lengths  similar  to  the  cross-shore  cuts  for  example,  might  result  in  the  identification  of 
more  of  these  areas  ) 
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Figure  56.      Cross-shore  Line  2«S,  Contributions  to  Density  Variance:     Contribution 

to  total  density  variance  (labeled  I )  by  temperature  (t),  salinity  (s),  and 
their  interaction  (t/s)  are  depicted  here  by  season  at  the  sin  lace  and  100 
m.  Units  of  density  variance  here  are  [kgm  xy. 
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Figure  57.  Alongshore  Line  6,  Contributions  to  Density  Variance:  Contribution 
to  total  density  variance  (labeled  T)  by  temperature  (1).  salinity  (s),  and 
their  interaction  (t/s)  are  depicted  here  by  season  at  the  surface  and  1<K) 
m.  Units  of  density  variance  here  arc  (kg  m  1)?. 
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Figure  58.      Alongshore  Line  3,  Contributions  to  Density  Variance:     Contribution 

to  total  density  variance  (labeled  T)  by  temperature  (t),  salinity  (s),  and 
then  interaction  (t/s)  are  depicted  here  by  season  at  the  surface  and  100 
m.  Units  of  density  variance  here  arc  (kg  nr3)2. 
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Figure  59.  Alongshore  Line  4,  Contributions  to  Density  Variance:  Contribution 
to  total  density  variance  (labeled  T)  by  temperature  (t),  salinity  (s),  and 
their  interaction  (t/s)  arc  depicted  here  by  season  at  the  surface  and  100 
ni.  Units  of  density  variance  here  are  (kg  m'3)2. 
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Figure  60.      Cross-shore  Line  27,  Contributions  to  Density  Variance:     Contribution 

to  total  density  variance  (labeled  T)  by  temperature  (t),  salinity  (s),  and 
their  interaction  (t  s)  are  depicted  here  by  season  at  the  surface  and  100 
m.  Units  of  density  variance  here  are  (kg  m  y)\ 
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Figure  61.      Cross-shore  Line  29,  Contributions  to  Density  Variance:     Contribution 

to  total  density  variance  (labeled  T)  by  temperature  (t),  salinity  (s),  and 
their  interaction  (t/s)  arc  depicted  here  by  season  at  the  surface  and  100 
m.  Units  of  density  variance  here  are  (kg  wr3)2. 


V.     SUMMARY  AND  CONCLUSIONS 

Salinity  has  a  profound  impact  on  the  CCS,  with  structure  that  is  complex  and 
variable.  In  the  region  of  the  Columbia  River's  outflow,  surface  salinity  may  vary  an 
entire  psu  in  100  km,  or  1  psu  over  the  course  of  a  year;  yet  in  this  region  the  depth  in- 
tegrated importance  of  the  salinity  gradient  is  less  important  than  anywhere  else  in  the 
CCS.  OIT  the  coast  of  Baja  California,  where  thermal  expansion  coefficients  are  high, 
saline  contraction  coefficients  low,  and  surface  salinity  gradients  have  minimal  impact 
on  the  density  field,  the  importance  of  salinity  variability  100  m  below  the  surface  is  of- 
ten equal  to  the  importance  of  temperature  variance.  The  mean  salinity  field  has  been 
shown  to  have  a  pronounced  effect  on  the  large-scale  geostrophic  circulation  of  the  CCS. 
On  the  other  hand,  there  seem  to  be  regions  on  the  order  of  3°  squares  where  the  as- 
sumption of  constant  salinity  might  not  cause  significant  errors  in  the  estimation  of 
density  gradients.  Relative  importance  of  temperature  and  salinity  often  changes  sea- 
sonally. Areas  where  mean  gradients  in  temperature  and  salinity  have  positive  corre- 
lation in  terms  of  their  effect  on  density  border  areas  of  negative  covariance  of  the  T-S 
gradients.  Despite  the  existence  of  correlation  in  mean  gradients,  there  seems  to  be  little 
correlation  between  simultaneously  observed  values  of  temperature  and  salinity. 

Based  on  these  observations,  it  is  reasonable  to  conclude  that  careful  work  con- 
cerning the  CCS  should  include  haline  description.  Descriptively  and  dynamically, 
salinity  is  essential  to  accurately  characterize  the  large-scale  structure  of  the  CCS.  Evi- 
dence has  been  developed  in  this  study  that  indicates  salinity  may  play  significant  roles 
in  many  smaller  scale  phenomenon  as  well.  There  is  some  indication  that  there  may  be 
areas  where  salinity  variation  may  be  neglected  in  process-oriented  models,  but  regions 
to  be  so  modeled  must  be  carefully  selected. 

A.     STERIC  ANOMALY 

A  comparison  of  the  constituents  of  the  steric  anomaly  due  to  salinity  and  temper- 
ature showed  that  the  gradients  in  <5r  and  Ss  oppose  generally,  major  exceptions  being 
on  the  surface  off  the  coast  of  Oregon  where  the  warm,  fresh  discharge  of  the  Columbia 
River  leads  to  large  values  of  o:  and  <55,  and  off  the  coast  of  northern  California  in 
summer,  when  upwelling  of  cold,  salty  water  leads  to  low  values  for  both  S7  and  Ss, 


In  areas  where  5T  is  highly  correlated  with  d,  and  has  similar  variance,  6T  might  be 
considered  an  adequate  estimator  of  3.  A  small,  shallow  region  with  these  characteristics 
was  identified  off"  the  coast  of  northern  Baja  California. 

Using  the  idea  of  RMSE  to  calculate  mean  differences  from  the  steric  anomaly,  S, 
of  Ss  and  <5r,  the  contributions  to  mean  gradient  in  d  associated  with  6S  and  ST  were 
compared.  On  the  large  scale,  <5S  is  quite  important  and  certainly  should  not  be  ignored. 
On  the  other  hand,  the  contribution  of  salinity  is  fairly  evenly  distributed  and  so,  as  the 
amount  of  area  considered  decreases,  the  relative  magnitude  of  error  in  gradient  resulting 
from  the  assumption  of  constant  salinity  also  decreases.  For  example,  assumption  of  a 
constant  value  of  salinity  for  the  entire  CCS  in  July  would  result  in  an  approximate 
overestimate  of  40  percent  in  the  alongshore  south  to  north  mean  6  gradient.  However, 
if  only  the  area  adjacent  to  the  California  coast  is  considered,  a  10  percent  error  results. 

B.  DYNAMIC  HEIGHT 

Calculating  dynamic  height,  <t>,  by  integrating  S  over  depth,  and  comparing  it  to 
Ojs,  ,  calculated  by  integrating  6T  ,  provided  insight  concerning  the  importance  of 
salinity  variation  to  the  geostrophic  flow  field.  At  all  depths,  the  mean  variability  of 
salinity  in  the  CCS  is  responsible  for  a  significant  equatorward  component  along  the 
coast  of  California  and  a  strong  offshore  component  adjacent  to  Baja.  Additionally 
mean  local  salinity  gradients  can  play  a  major  role  in  the  structure  of  small  scale 
features. 

C.  VARIANCE  IN  THE  DENSITY  FIELD 

By  calculating  temperature  and  salinity  variances  and  covariances  and  differentiating 
the  equation  of  state,  contributions  to  local  variation  in  density  due  to  salinity,  temper- 
ature, and  their  interaction  were  computed  directly.  Temperature  variation  provided  the 
major  contribution  year-round  on  the  surface.  Salinity's  importance  generally  increased 
with  depth  so  that  at  500  m,  salinity  variance  might  actually  be  most  significant.  Sub- 
stantial variance  in  the  density  field  existed  even  at  500  m  in  locations  that  correspond 
to  the  CC  and  CU,  indicating  dynamic  activity  was  occurring  at  500  m,  and  perhaps 
that,  as  Lynn  and  Simpson  (1989)  have  suggested,  the  CCS  can  extend  deeper  than 
500  m.  and  that  at  deeper  depths  it  may  be  traceable  by  salinity  variance. 

By  applying  this  technique  to  the  mean  temperature  and  salinity  field,  the  relative 
contribution  to  longshore  and  cross-shore  gradient  was  determined.  Nearshore.  or 
elsewhere  in  the  vicinity  of  river  drainage  (Figure  52).  variance  is  generally  high  and 
saline  contribution  considerable.    There  is  a  tendency  for  hich  variance  to  occur  in  the 


autumn  (defined  here  as  October  through  December,  Figure  57)  due  apparently  to  sur- 
face cooling.  Offshore  this  cooling  and  associated  mixing  tend  to  decouple  the  surface 
layer  and  allow  greater  density  variation  at  deeper  depths,  at  least  100  m. 


VI.      RECOMMENDATIONS 

A.  PROFILE 

The  information  developed  in  this  research  should  be  displayed  in  vertical  sections. 
The  organizational  techniques  employed  here  would  allow  alongshore,  cross-shore,  or 
diagonal  lines  to  be  displayed  easily  in  vertical  cross  section.  Study  of  these  sections 
would  allow  fuller  appreciation  of  the  three  dimensional  interaction  of  temperature  and 
salinity  in  the  CCS.  Perhaps  sections  of  the  contribution  of  salinity  to  density  variance 
could  be  used  with  the  ideas  of  Lynn  and  Simpson  (1989)  concerning  "spiciness"  to  gain 
a  fuller  understanding  of  the  water  mass  structure  of  the  CCS. 

B.  FINER  RESOLUTION 

This  analysis  has  been  carried  out  at  a  1°  spatial  by  I  month  temporal  scale,  and  so 
allows  only  large-scale  resolution.  A  much  finer  scale  approach,  perhaps  1,2°  by  2 
weeks,  might  be  possible  in  certain  regions.  Using  similar  processing  techniques  and 
finer  resolution,  if  even-  depth  level  were  analyzed,  rather  than  the  few  depths  empha- 
sized here,  the  T-S  interplay  associated  with  small  scale  phenomena  such  as  coastal 
upweiling  might  be  characterized. 

C.  MODELING  CONSIDERATIONS 

Speculation  has  been  offered  here  concerning  the  validity  of  the  assumption  of 
constant  salinity  for  the  purpose  of  dynamic  modeling.  Ideally,  the  mean  salinity  data 
processed  for  this  study  should  be  incorporated  into  a  CCS  model  run  and  results  com- 
pared to  similar  runs  in  which  salinity  was  held  constant.  On  the  other  hand,  indications 
of  the  effects  of  salinity  might  be  achieved  with  less  effort  by  using  a  temperature-driven 
linearization  of  density  that  simulates  saline  variability.  For  example,  increasing  the 
magnitude  of  the  thermal  expansion  coefficient,  a,  in  an  upweiling  process  model  may 
capture  the  effect  of  a  positive  T-S  covariance  on  density.  Entering  different  values  of 
a  for  different  regions  and  depths  (higher  values  where  covariance  is  positive,  lower  val- 
ues where  negative)  would  simulate  salinity's  variation  over  regions  where  salinity  or 
temperature  effects  change  sign.  Use  of  a  constant  value  of  a  over  a  range  of  temper- 
atures, in  addition  to  ignoring  saline  contribution,  can  superimpose  a  false  smoothing 
tendency,  resulting  in  an  underestimation  of  the  magnitude  of  positive  Ap  in  cold  water, 
and   negative   Ap   in   warm   water.      Forcing   a   to   change   more    or   less   uniformly 


meridionally  but  at  a  higher  rate  than  warranted  by  the  mean  temperature  gradient 
would  allow  study  of  the  combined  effects  of  meridional  saline  variation  and  smoothing. 
Because  salinity  measurements  are  more  expensive  to  make  and  more  difficult  to 
obtain  than  temperature  observations,  it  may  be  desirable  to  estimate  salinity  from  a 
known  temperature-salinity  (T-S)  correlation  (Stommel,  1947).  Emery  (1975)  has  shown 
that  the  accuracy  of  dynamic  height  calculations  made  from  such  derived  salinities  can, 
at  least  in  some  cases,  be  comparable  to  the  accuracy  of  determinations  made  from  ob- 
served salinities.  In  these  cases,  the  salinity  field  in  general  can  be  accurately  inferred 
from  expendable  bathythermograph  observations  alone,  as  Rienecker  et  al.  (1985)  have 
attempted. 

The  validity  of  such  salinity  inferences  depends  completely  on  the  strength  of  the 
T-S  correlation,  or  what  Stommel  (19-47)  called  "range  of  uncertainty",  i.e.,  the 
"tightness"  of  a  particular  T-S  curve.  This  uncertainty  in  differential  salinity  can  be  re- 
lated to: 


N./'(l-r*)x(AS)',  (17) 

where  AS  is  determined  from  the  T-S  curve,  rts  is  the  coefficient  of  determination  for 
temperature  and  salinity,  and  (1  -  r-s)  is  the  fraction  of  salinity  variance  that  cannot  be 
explained  by  temperature  variance.  Clearly  the  uncertainty  increases  as  the  strength  of 
the  correlation  (and  magnitude  of  r*)  decreases.  Therefore,  it  may  be  reasonable  to  ex- 
pect the  method  of  salinity  inference  from  T-S  relationship  to  be  most  accurate  in  re- 
gions of  highest  T-S  correlations,  and  least  accurate  in  regions  where  the  T-S  correlation 
is  weak.  Blumberg  et  al.  (19S-4)  concluded  that  development  of  T-S  curves  for  the  CCS 
would  not  be  cost  effective  for  modeling  eilorts  due  to  combinations  of  relatively  poor 
T-S  correlations  inherent  in  near  surface  waters  (less  than  200  m),  and  lack  of  significant 
numbers  of  temperature  observations  in  deeper  areas  where  the  idea  might  otherwise 
work. 

Preliminary  investigations  have  been  made  to  determine  seasonal  values  of  the  T-S 
correlation  (i.e..  coefficients  of  determination)  at  10,  100.  200  and  500  m.  Even  at 
500  m,  areas  where  temperature  variance  has  as  much  as  a  50  percent  linear  relationship 
with  salinity  are  rare  (figure  62),  and  those  areas  are  not  constant  from  month  to 
month.  This  result  strongly  supports  the  statements  of  Blumberg  et  al.  ( 1 9S4).  A 
modeler  might  be  served  as  well  by  using  salinity  climatology  as  by  attempting  to  pre- 
scribe   the      saliniu       field      based      on   its   relationship   with    temperature.  An 


observationalist  should  probably  take  full  CTD  casts  to  determine  the  salinity  in  the 
CCS.    Calculations  could  be  performed  to  confirm  these  statements. 

The  geographic  distribution  of  salinity  in  the  CCS  may  be  sufficiently  linear,  at  least 
regionally,  to  allow  worthwhile  use  of  linear  regression  coefficients  on  latitude,  longitude 
and  or  depth.  Sums  required  by  least  squares  techniques  to  compute  statistically  sig- 
nificant coefficients  of  such  a  regression  equation  are  easily  obtainable  through  tech- 
niques used  in  this  study.  Model  runs  using  the  resulting  salinity  parameterization 
would  be  instructive,  as  would  the  shorthand  description  of  the  salinity  field  represented 
by  the  parameterization  itself.  If  results  of  any  of  these  model  runs  differ  significantly 
from  corresponding  constant  salinity  runs,  decisions  concerning  the  best  method  to 
include  salinity  could  be  made. 

D.     SOUND  SPEED  CALCULATIONS 

Salinity  is  a  variable  in  the  equation  for  sound  speed  in  seawater  but  is  generally 
estimated  by  a  climatological  mean  when  tactical  acoustic  range  predictions  are  com- 
puted. The  considerable  database  now  available  in  the  CCS  should  allow  a  numerical 
evaluation  of  the  validity  of  such  a  parameterization. 
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Figure  62.  Coefficient  of  Determination  (T-S),  Slimmer,  500  m:  r*  (the  algebraic 
sign  is  obtained  from  covariance,  and  "99"  indicates  a  bad  data  proc- 
essing point). 


APPENDIX  A.     NUMBER  OF  OBSERVATIONS 

This  appendix  provides  plots  of  the  number  of  observations  associated  with  each 
data  block  (for  those  data  blocks  having  at  least  ten  observations)  for  the  months 
January,  April.  July,  and  October  and  the  depths  10  and  500  m. 
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Figure  63.      Number  of  Observations  at  10  in,  January. 
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Figure  64.      Number  of  Observations  at  10  m,  April. 
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Figure  65.      Number  of  Observations  al  10  in,  July. 
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Figure  66.      Number  of  Observations  at  10  m,  October. 
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Figure  67.      Number  of  Observations  at  500  in,  J  a  mini  v. 
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Figure  68.      Number  of  Observations  at  500  in,  April. 
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Figure  69.      Number  of  Observations  nl  500  in,  July. 
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Figure  70.       Number  of  Observations  at  500  in,  October. 
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APPENDIX  B.     SMALL  SCALE  EFFECTS  OF  SALINITY  ON  DYNAMIC 
HEIGHT 

In  addition  to  introducing  the  general  large-scale  biases  to  the  estimated 
geostrophic  flow  field  discussed  in  section  IVE,  the  OiSrp  estimator  may  misrepresent 
some  smaller  scale  features.  The  December  mean  O  map  shows  a  strong  closed  cyclonic 
circulation  near  30°  N,  121°  W.  The  December  mean  <£>iSrp  field  also  shows  a  strong 
cyclonic  feature  at  that  point,  indicating  a  thermal  front  as  the  probable  cause  of  the 
enhanced  flow.  However,  the  difference  map  for  the  period  also  places  a  cyclonic  cir- 
culation at  that  location,  showing  the  front  to  be  haline  as  well  as  thermal. 

In  April  a  closed  cyclone  again  appears  on  the  O1  map.  this  time  near  32°  N, 
123°  \V.  Here  the  d>35;.,  map  shows  only  weak  cyclonic  turning  in  basically  onshore  flow. 
The  corresponding  cp  map  indicates  strong  horizontal  cyclonic  shear  in  the  region.  In 
this  case,  changes  in  salinity  over  a  relatively  small  distance  induce  considerable  shear 
in  the  mean  field.  This  feature  would  be  missed  entirely  by  an  estimate  based  on  tem- 
perature variation  alone. 

In  May  a  similar  cyclonic  anomaly  appears  on  the  Ox  map  near  36°  N,  127°  W. 
The  <I)  field  shows  a  strong  ridge  at  37°  30v  N  and  a  strong  trough  at  36°  X,  while  the 
(f  ,,.  map  shows  a  trough  at  37.5°N  and  a  ridge  at  36.5°\.  Admittedly  the  scale  of  this 
mean  feature  is  near  the  resolvable  limit,  and  some  critical  mean  values  of  O  have  been 
determined  from  relatively  few  casts.  However,  if  the  feature  is  correctly  represented,  it 
demonstrates  the  importance  of  salinity  as  a  modifier  of  the  effects  of  temperature.  It 
is  interesting  that  O  (200  500)  and  <I>35v  (200,  500)  contours,  by  contrast,  are  in  phase  in 
this  area.  Analysis  of  the  May  salinity  distribution  over  the  area  of  interest  shows  sharp 
salinity  gradients  at  the  surface  and  less  pronounced  gradients  at  depth.  Consequently, 
the  inferred  flow  fields  are  consistent. 
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